Bose-Einstein condensation and Superfluidity of magnetoexcitons in Graphene 
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We propose experiments to observe Bose-Einstein condensation (BEC) and superfluidity of quasi- 
two-dimensional (2D) spatially indirect magnetoexcitons in bilayer graphene. The magnetic field 
B is assumed strong. The energy spectrum of collective excitations, the sound spectrum as well as 
the effective magnetic mass of magnetoexcitons are presented in the strong magnetic field regime. 
The superfluid density ns and the temperature of the Kosterlitz-Thouless phase transition T c are 
shown to be increasing functions of the excitonic density n but decreasing functions of B and the 
interlayer separation D. Numerical results are presented from these calculations. 

PACS numbers: 71.35.Ji, 71.35.Lk, 71.35.-y 



Indirect excitons in coupled quantum wells (CQWs) 
in the presence or absence of a magnetic field B have 
been the subject of recent experimental investigations 
EL B S H- These systems are of particular interest 
because of the possibility of Bose-Einstein condensation 
(BEC) and the superfluidity of indirect excitons formed 
from electron-hole pairs. These may result in persistent 
electrical currents in each Q W or coherent optical prop- 
erties and Josephson junction phenomena @j @, 0i fi Q- 
In high magnetic fields, two-dimensional (2D) excitons 
survive in a substantially wider temperature range, as 
the exciton bindin g en ergies increase with magnetic field 

In this Letter we propose a new physical realization 
of magnetoexcitonic BEC and superfluidity in bilayer 
graphene with spatially separated electrons and holes 
in high magnetic field. Recent technological advances 
have allowed the production of graphene, which is a 2D 
honeycomb lattice of carbon atoms that form the basic 
planar structure in graphite E3, EH Graphene has been 
attracting a great deal of experimental and theoretical 
attention because of unusual properties in its bandstruc- 
ture El) US HH H2] • It is a gapless semiconductor with 
massless electrons and holes which have been described 
as Dirac-fermions [23j . Since there is no gap between 
the conduction and valence bands in graphene without 
magnetic field, the screening effects result in the absence 
of excitons in graphene in the absence of magnetic field. 
A strong magnetic field produces a gap since the energy 
spectrum becomes discrete formed by Landau levels. The 
gap reduces screening and leads to the formation of mag- 
netoexcitons. 

We consider two parallel graphene layers separated by 
an insulating slab of Si02 . The electrons in one layer and 
holes in the other can be controlled as in the experiments 
with CQWsEL H S H by laser pumping (far infrared in 
graphene). The spatial separation of electrons and holes 
in different graphene layers can be achieved by apply- 
ing an external electric field. Furthermore, the spatially 
separated electrons and holes can be created by varying 
the chemical potential by using a bias voltage between 



two graphene layers or between two gates located near 
the corresponding graphene sheets. Indirect magnetoex- 
citons are bound states of spatially separated electrons 
and holes in an external magnetic field. The ratio of the 
external voltage V ex t to the interlayer separation D re- 
quired to create spatially separated electrons and holes 
in graphene layers with the 2D density n — 10 1:l cto -2 is 
given by V ext /D = ^enD/t b = 4.021 x 10 4 V/cm. Here, 
— e is the electron charge and = 4.5 is the dielectric 
constant of Si02- Since the critical electric field E cr of 
the dielectric breakdown for Si02 is E cr as 10 6 V/cm, we 
conclude that the external electric field for the spatially 
separated electrons and holes is less than the critical elec- 
tric field for dielectric breakdown in Si02- 

A conserved quantity for an isolated electron-hole pair 
for graphene in magnetic field B is the exciton magnetic 
momentum P defined as 



P = -iftV a - ihV h + -(A e - A h ) - -[B x (r e - r h )\ (1) 

c c 

for the Dirac equation[24| as for the Schrodinger equation 
EH E2, HH- Here, r e and are 2D coordinate vectors 
of an electron and hole, respectively, A e and A^ are the 
corresponding vector potential of an electron and hole. 
The cylindrical gauge for vector potential is used with 
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Neglecting transitions between Landau levels for high 
magnetic fields, we employ first-order perturbation the- 



ory to the Coulomb attraction V(r) 



2 /(eWr 2 + D*) 



between an electron and hole. Here, r = r e rv We cal- 
culate the magnetoexciton energy using the expectation 
value for an electron in Landau level I and a hole in level 
1. We have 
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where \h, m, P) is an eigenfunction of the non-relativistic 
Hamiltonian of a non-interacting electron-hole pair de- 
fined in [ToL Il4j | . For small magnetic momentum satisfy- 
ing P <^hJjB and P <C HD/r B , we obtain the following 
relations [14{ 



((nmP\V(r)\nmP)) = E\ 



hm 
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2M ilm (B,D) ' 

Substituting these results into Eq. we obtain the 

binding energy £ B {D) and the effective magnetic mass 
m B (D) of a magnetoexciton with spatially separated 

electron and hole in bilayer graphene as £ b \d) = 
£&\B,D) + 2£$(B,D) + E$(B,D), and m B \D) = 
M^(B,D) + MoiO^-D)" 1 + M^(B,D), where con- 
stants £${B,D), E$(B,D), 8§\B,D), M 00 (B,D), 
Mqx(B,D") and M 10 (B,D) depending on magnetic field 
B and the interlayer separation D in detail are given in 
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Here, tb — yJftc/{eH) is the magnetic length; c 
is the speed of light, £ = e 2 /ebi~B M — 

2 3 / 2 S 2 ef,/(v / 7re 2 rs) and erfc(z) is the complementary er- 
ror function [T3| ■ The radius of a magnetoexciton in the 
lowest Landau level is given by n i i(B) = 4re. 

For large interlayer separation D ^> rs, the asymp- 
totic values of the binding energy £ { b\d) and the effec- 
tive magnetic mass ms (D) axeSg (B,D) = —4e 2 /(e b D), 
and m B {D) = e b D 3 B 2 /(4c 2 ). When D < r B , these 
quantities denoted by £q and Mo are presented above. 
We can see that the effective magnetic mass of an indirect 
magnetoexciton is approximately four times smaller than 
in CQWs at the same D, eb and B [14] . The magnetoexci- 
tonic energy is approximately four times larger in bilayer 
graphene than in CQWs. Measuring energies relative to 
the binding energy of a magnetoexciton, the dispersion 
relation £fc(P) of a magnetoexciton is quadratic at small 
magnetic momentum, i.e., P <C h/rs and P <C hD/r 2 B . 
We have £fe(P) = P 2 /(2rriBk), where rriBk is the effective 
magnetic mass, dependent on B, D and the magnetoex- 
citonic quantum numbers k = {n + ,n_} for an electron 
in Landau level n+ and a hole in level n_. Indirect mag- 
netoexcitons, cither in the ground state or an excited 
state, have electrical dipole moments. We treat these 
excitons as interacting parallel electric dipoles. This is 
valid when D is larger than the mean separation (r) be- 
tween an electron and hole parallel to the graphene lay- 
ers. We take into account that at high magnetic fields 



(r) w Pr 2 B /h with (r) perpendicular to P. Typical val- 
ues of magnetic momenta are given by P ~ fry/n, where 
n is the 2D density of magnetoexcitons for a parabolic 
dispersion relation. Consequently, D (r) is valid when 
D ^S> y/rirg . Since electrons on a graphene lattice can be 
in two valleys, there are four types of excitons in bilayer 
graphene. Due to the fact that all these types of excitons 
have identical envelope wave functions and energies [24|. 
we consider below only excitons in one valley. Also, we 
use rio = n/4 as the density of excitons in one layer, 
with n denoting the total density of excitons. For large 
electron- hole separation D r Bl transitions between 
Landau levels due to the Coulomb electron-hole attrac- 
tion can be neglected, if the following condition is valid, 
i.e., Eb = 4e 2 /(e;,Z?) <C hvp/r B . This corresponds to 
high magnetic field B, large interlayer separation D and 
large dielectric constant of the insulator layer between 
the graphene layers. In this notation, vp — \/3at/(2h) is 
the Fermi velocity of electrons. Also, a — 2.566A is a lat- 
tice constant, t w 2.7leV is the overlap integral between 
nearest carbon atoms [26| . 

The distinction between excitons and bosons is due to 
exchange effects @. These effects for excitons with spa- 
tially separated electrons and holes in a dilute system 
satisfying na 2 (B, D) < 1 are suppressed due to the neg- 
ligible overlap of the wave functions of two excitons as a 
result of the potential barrier, associated with the dipole- 
dipole repulsion [9( • Two indirect excitons in a dilute sys- 
tem interact via U(R) — e 2 D 2 / (e b R 3 ), where R is the dis- 
tance between exciton dipoles along the graphene layers. 
In high magnetic fields, the small parameter mentioned 
above has the form exp[— 2h~ 1 (m,Bk) 1 ^ 2 eDa~ 1 / 2 (B, D)\. 
So at T — 0, the dilute gas of magnetoexcitons, which is 
a boson system, form a Bose condensate [27j ■ Therefore, 
the system of indirect magnetoexcitons can be treated by 
the diagrammatic technique for a boson system. For the 
dilute 2D magnetoexciton system with na 2 (B,D) <C 1, 
the sum of ladder diagrams is adequate. For the low- 
est Landau level, we denote £n(P) = e(P)- Using 
the orthonormality of the four-component wave functions 
of the relative coordinate for a non-interacting pair of 
an electron in Landau level n+ and a hole in level n_ 
(^ > n + ,n_(P = 0,r))[24j we obtain an approximate equa- 
tion for the vertex T in strong magnetic fields. Due 
to the orthonormality of the four-component wave func- 
tions $ ra+iTl _ (0, r) the projection of the equation for the 
vertex in the ladder approximation for a dilute system 
onto the lowest Landau level results in the scalar inte- 
gral equation which does not reflect the spinor nature of 
the four-component magnetoexcitonic wave functions in 
graphene. In high magnetic field, one can ignore transi- 
tions between Landau levels and consider only the lowest 
Landau level states n + = n_ = 1. Since typically, the 
value of r is rs, and P <C fr/rs in this approximation, 
the equation for the vertex in the magnetic momentum 
representation P for the lowest Landau level has the same 
form (compare with ||) as for a 2D boson system in the 
absence of magnetic field, but with the magnetoexciton 
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magnetic mass raj (which depends on B and D) instead 
of the exciton mass (M = m e + ms) and magnetic mo- 
mentum instead of inertial momentum: 
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n Q T = n o r(0, 0;0) , 
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where P = {P, f2}, and fj, is the chemical potential of the 
system. Equation ([5]) is valid at parameter values which 
satisfy the condition for validity of the perturbation the- 
ory applied to the calculation of the magnetoexcitonic 
binding energy. The specific feature of a 2D Bose sys- 
tem is connected with a logarithmic divergence in the 
2D scattering amplitude at zero energy 0, |{|. A sim- 
ple analytical solution of Eq. (J5]) for the chemical poten- 
tial can be obtained if KmBe 2 D 2 /(h 3 e) <C 1. In strong 
magnetic fields at J) > the exciton magnetic mass 
is defined as tub ~ h 2 eD 3 /(4e 2 r%). So the inequality 
nm B e 2 D 2 /(h 3 e) < 1 is valid if D < (r 4 /n 1 / 2 ) 1 / 5 . Con- 
sequently, the chemical potential /i is obtained as 



/' 



2m B 8m B log[ft 4 e 2 /(27rnm 2 3 e 4 D 4 )]' 
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At small magnetic momentum, the solution of Eq.(f5|) 
corresponds to the sound spectrum of collective exci- 
tations e(P) — c s P. Here, the sound velocity c s = 
y/Tn/ (4ms) = y /V m B; where /i is guven by Eq. ([6]). 
Since magnetoexcitons have a sound spectrum of col- 
lective excitations at small magnetic momentum P due 
to the dipole-dipole repulsion, the magnetoexcitonic su- 
perfluidity is possible at low temperature T in bilayer 
graphene. This is so since the sound spectrum satisfies 
the Landau criterium of superfluidity [271 ] . 

It can be shown that when D = 0, the interac- 
tion between two magnetoexcitons in the lowest Lan- 
dau level can be neglected in strong magnetic field [l(|. 
The magnetoexcitons constructed by spatially separated 
electrons and holes in bilayer graphene at large inter- 
layer separations D 3> re form a weakly interacting 
2D non-ideal Bose gas with a dipole-dipole repulsion. 
Thus, the phase transition from the normal to super- 
fluid phase is the Kosterlitz-Thouless transition [28| . The 
temperature of this phase transition T c to the superfluid 
state in a 2D magnetoexciton system is determined from 
T c = Trh 2 ns{T c )/ (2kBrriB), where ns(T) is the superfluid 
density of the magnetoexciton system, as a function of 
T, B, D and ks is Boltzmann's constant. The func- 
tion ns{T) can be obtained from ns = n/ '4 — n n , with 
n the total density and njv the normal component den- 
sity. To calculate the superfluid component density, we 
find the total quasiparticle current in a reference frame 



in which the superfluid component is at rest. We deter- 
mine the normal component density by the usual pro- 
cedure [27|. Suppose that the magnetoexciton system 
moves with a velocity u. At nonzero temperatures dis- 
sipating quasiparticles will appear in this system. Since 
their density is small at low temperatures, one can as- 
sume that the gas of quasiparticles is an ideal Bose gas. 
To calculate the superfluid component density, we find 
the total current of quasiparticles in a frame of reference 
in which the superfluid component is at rest. We de- 
note (Pk\...\Pk) as (•••). Using the Feynman theorem 
for isolated magnetoexcitons, we obtain the velocity [25[ 

v = (v) = (dH/dP) = de k (P)/dP = P/m Bk . We 



obtain the mean total current of 2D magnetoexcitons 
in the coordinate system, moving with a velocity u as 
(J) = (P) /tub- Expanding the integrand to first order 
by Pu/(fcfiT), we have 
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where f[e(P)} = (exp[e{P)/(k B T)} - 1) 1 is the Bose- 
Einstein distribution function and £(z) is the Riemann 
zeta function (£(3) ~ 1.202). Then we define the nor- 
mal component density is [2 71] (J) = n n u. Applying 
Eq. ([7]), we obtain the expression for the normal den- 
sity tin- As a result, we have for the superfluid density: 
ns{T) = n/A — jin(T). It follows that the expression for 
the superfluid density ns in strong magnetic field for the 
proposed magnetoexciton system differs from the analo- 
gous expression in the absence of magnetic field in semi- 
conductor CQWs (compare with Refs. [9() by replacing 
the total exciton mass M = m e + tub with the magne- 
toexciton magnetic mass ms- 

In a 2D system, superfluidity of magnetoexcitons ap- 
pears below the Kosterlitz-Thouless transition tempera- 
ture, where only coupled vortices are present [281 ]. Em- 
ploying ns(T) for the superfluid component, we obtain 
an equation for the Kosterlitz-Thouless transition tem- 
perature T c with solution 
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Here, is an auxiliary quantity, equal to the tem- 
perature at which the superfluid density vanishes in 
the mean- field approximation, i.e., ns(T®) = 0, = 

kg 1 (27rft 2 nc 4 TO B /(12C(3))) 1/3 . The temperature T c ° = 
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FIG. 1: (Color on line) Dependence of Kosterlitz-Thouless 
transition temperature T c — T C (B) (in units K) versus mag- 
netic field for bilayer graphene separated by Si02- with 
£b = 4.5. The magnetoexciton density n = 4 x lO^cm -2 . 
Different interlayer separations D are chosen: D — 30nm 
(solid curve), D = 28nm (dotted curve), D = 27nm (dashed 
curve) . 

T®(B,D) may be used to estimate the crossover region 
where local superfluid density appearers for magnetoex- 
citons on a scale smaller or of the order of the mean in- 
tervortex separation in the system. The local superfluid 
density can manifest itself in local optical or transport 
properties. The dependence of T c on B and D is repre- 
sented in Fig. [1] According to Eq. © , the temperature 
T c for the onset of superfluidity due to the Kosterlitz- 
Thouless transition at a fixed magnetoexciton density de- 
creases as a function of magnetic field B and interlayer 
separation D. This is due to the increased tub as a func- 



tions of B and D. The T c decreases as B 1 / 2 at D <C tb 
or as B~ 2 when D ^S> Tu- 
la conclusion, we have studied BEC and superfluid- 
ity of magnetoexcitons in two graphene layers with ap- 
plied external voltage in perpendicular magnetic field. 
The superfluid density ns{T) and the temperature of 
the Kosterlitz-Thouless phase transition to the super- 
fluid state have been calculated. We have shown that at 
fixed exciton density n the Kosterlitz-Thouless tempera- 
ture T c for the onset of superfluidity of magnetoexcitons 
decreases as a function of magnetic field like B^ 1 / 2 at 
D <tb and as B~ 2 when D tb- We have shown that 
T c increases when the density n increases and decreases 
when the magnetic field B and the interlayer separation 
increase. The dependence of T c on B and D is presented 
in Fig. [TJ We also note that the superfluidity of indi- 
rect magnetoexcitons in strong perpendicular magnetic 
field in bilayer graphene is very interesting, because the 
magnetoexcitons in graphene are found to be more stable 
than in CQWs. Namely, the binding energy of magne- 
toexcitons in graphene is four times greater than that in 
CQWs with the same D, ej, and B. We consider only the 
collective properties of excitons with electrons and holes 
from the same valley. We note that there is no crossover 
between Bose condensates of different types of excitons. 
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